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ABSTRACT Silene latifolia is a dioecious plant with heteronnorphic XY sex chronnosonnes. Previous studies 
of sex chronnosonne-linked genes have suggested a gradual divergence between the X-linked and the Y-linked 
genes in proportion to the distance from the pseudoautosomal region. However, such a comparison has yet to be 
made for the noncoding regions. To better characterize the nonrecombining region of the X and Y chromosomes, 
we sequenced bacterial artificial chromosome clones containing the sex chromosome-linked paralogs SIAP3X 
and SIAP3Y, including 115 kb and 73 kb of sequences, respectively, flanking these genes. The synonymous 
nucleotide divergence between SIAP3X and SIAP3Y indicated that recombination stopped approximately 3.4 
million years ago. Sequence homology analysis revealed the presence of six long terminal repeat retrotransposon- 
like elements. Using the nucleotide divergence calculated between left and right long terminal repeat sequences, 
insertion dates were estimated to be 0.083-1.6 million years ago, implying that all elements detected were 
inserted after recombination stopped. A reciprocal sequence homology search facilitated the identification of 
four homologous noncoding DNA regions between the X and Y chromosomes, spanning 6.7% and 10.6% of the 
X chromosome-derived and Y chromosome-derived sequences, respectively, investigated. Genomic Southern 
blotting and fluorescence in situ hybridization showed that the noncoding DNA flanking SIAP3XA^ has homology 
to many regions throughout the genome, regardless of whether they were homologous between the X and Y 
chromosomes. This finding suggests that most noncoding DNA regions rapidly lose their counterparts because of 
the introduction of transposable elements and indels (insertion-deletions) after recombination has stopped. 
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Many plants, including the model plant Arahidopsis thaliana, are bi- 
sexual. However, approximately 6% of angiosperms are dioecious, 
meaning that male plants produce only male flowers and female 
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plants produce only female flowers (Renner and Ricklefs 1995). Some 
dioecious plants have sex chromosomes with varying sizes and forms. 
Carica papaya (60 Mb) (Ming et al 2008) and Asparagus officinalis 
(65 Mb) (Telgmann-Rauber et al 2007), for example, have homomor- 
phic sex chromosomes. In contrast, Marchantia polymorpha (X: 20 
Mb; Y: 10 Mb) (Yamato et al 2007) and Rumex acetosa (X: 740 Mb; 
Yi: 520 Mb; Y2: 410 Mb) (Mosiolek et al 2005) have heteromorphic 
sex chromosomes. In Silene latifolia, sex is determined by heteromor- 
phic sex chromosomes, with a pair of X chromosomes in females and 
both an X chromosome and a Y chromosome in males. The Y chro- 
mosome is estimated to be 570 Mb (Siroky et al 2001; Liu et al 2004), 
representing approximately 9% of the total genome (Matsunaga et al 
1994), and is 1.4-times larger than the X chromosome (Matsunaga 
et al 1994). Thus, the X and Y chromosomes of S. latifolia provide 
a good opportunity to study the evolutionary history of heteromor- 
phic sex chromosomes during plant evolution. 
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Sex chromosomes have evolved independently in many plant 
groups (Charlesworth 2002). Sex chromosomes are thought to have 
originated from a pair of autosomes with two sexually antagonistic 
mutations in plants and in animals. It is also assumed that the chro- 
mosomal region harboring these mutations was subject to selection- 
driven suppression of recombination (Charlesworth 2013), which may 
have been facilitated by chromosomal inversions (Lemaitre et al. 2009; 
Wang et al 2012) and translocations (Charlesworth and Charlesworth 
1980). In mammals, such nonrecombining genomic regions are 
known to have expanded over evolutionary time (Iwase et al. 2003). 

One way to estimate the age of sex chromosomes is to study X-Y 
divergence (Charlesworth 2013). Estimating the date of sex chromo- 
some emergence is possible when using the sex-determining gene 
itself and the use of silent site divergence of an X-Y gene pair allows 
for estimation of the time when recombination stopped between the 
two chromosomal regions in which the gene pair resides. Using sev- 
eral X-Y gene pairs, dates associated with recombination suppression 
have been estimated in species of mammals (Lahn and Page 1999), 
birds (Lawson-Handley et al. 2004; Nam and Ellegren 2008), and 
plants [S. latifolia (Bergero et al. 2007) and C. papaya (Wang et al. 
2012)], revealing differences in the estimated dates based on the chro- 
mosomal regions in which the pairs were located. Chromosomal 
regions exhibiting different silent site sequence divergence estimates 
are referred to as "evolutionary strata" (Lahn and Page 1999). The 
oldest stratum in S. latifolia is thought to have appeared 5-10 million 
years ago (MYA) (Nicolas et al. 2005). By mapping eight X-linked 
genes, it has been reported that silent site divergence between X-Y 
gene pairs increased in proportion to the distance of the gene pair 
from the pseudoautosomal region (PAR), and it has been suggested 
that recombination between the X and Y chromosomes stopped in 
progressive steps that formed two evolutionary strata (Bergero et al. 
2007). Bergero et al. (2013) showed another stratum that was formed 
by additions of genome regions. However, previous studies focused on 
the comparison of only X-Y gene pairs or the accumulation of re- 
petitive sequences at the chromosomal level (Kejnovsky et al. 2009). 
Thus, the nonrecombining regions that comprise each evolutionary 
stratum remain uncharacterized. 

The S. latifolia MADS box gene, S/AP3, which exhibits significant 
similarity to A. thaliana APETALAS, was first isolated by Matsunaga 
et al. (2003) and was found to occur in two copies, SIAP3Y and 
SIAP3A. SIAP3Y hdis only been detected in male plants when genomic 
PCR is performed using S/AP3Y- specific primers and is thus known to 
be located on the Y chromosome. In contrast, SIAP3A was thought to 
be located on the autosome, because this gene was first detected by 
genomic PCR using S/AP3A-specific primers and flow-sorted X chro- 
mosomes and autosomes (Matsunaga et al. 2003). It was also later 
amplified by Cegan et al. (2010) by using microdissected X chromo- 
somes, but not microdissected autosomes; the gene has since been 
renamed SIAP3X. Reaching the same conclusion, Nishiyama et al. 
(2010) conducted a segregation analysis of SIAP3X in which different 
lines of S. latifolia were crossed and the segregation of the male 
parent-derived sequence was examined. The initial misunderstanding 
of SIAP3X mainly results from the use of a localization test of flow- 
sorted chromosomes derived from cultured root cells in which the 
translocation including SIAP3X would have occurred. Cegan et al. 
(2010) compared promoter sequences of SIAP3X and SIAP3Y and 
revealed that a specific sequence was inserted in the S/AP3Y promoter. 
This inserted sequence shows some homology to the promoter 
sequence of MROSl, which is specifically expressed in male plants 
(Matsunaga et al. 1996, 1997). The divergence of noncoding regions 
is possibly a driving force behind the degeneration of X and Y 



chromosomes, resulting in differential expression of X-linked and 
Y-linked genes. 

Based on the synonymous nucleotide divergence calculated in the 
coding regions of SIAP3X and SlAP3Yy we estimated that recombina- 
tion between the two genes stopped approximately 3.4 MYA. Com- 
parison of intron sequences within the pair of genes revealed that 
homologous regions have been segmentalized by insertions and dele- 
tions (indels) and contain repetitive sequences. Additionally, we com- 
pared sequences surrounding the gene pair; predictably, the transposons 
classes observed in these regions differed between the X and Y chro- 
mosomes. We estimated insertion times of predicted long terminal re- 
peat (LTR) retroelements based on the nucleotide divergence between 
left and right LTRs, which implied that they had likely been inserted in 
these regions after recombination between SIAP3X and S/AP3 7 stopped. 
These data provide the first evidence related to nucleotide resolution of 
sequence divergence between a pair of noncoding DNA regions on the 
sex chromosomes of S. latifolia after recombination stopped. 

MATERIALS AND METHODS 
Plant materials 

An inbred S. latifolia line, the K line (Kazama et al. 2003), was used 
for all experiments. A second line, the B line, was used for linkage 
mapping of the X chromosome. The original plants of K and B lines 
were provided by the University of Oslo Botanical Garden, Norway (K 
line), and the Royal Botanic Gardens, Kew, United Kingdom (B line). 
Fi progeny of the two phyletic lines were obtained by crossing a male 
plant of the K line with a female plant of the B line. F2 progeny were 
generated by crossing an Fi male with an Fi female. Plants were 
grown in pots in a regulated chamber at 23° under a 16:8-hr light/ 
dark cycle. Young leaves were used for genomic DNA isolation. 

Bacterial artificial chromosome sequencing 

The clone 2b5E (renamed as 13dllEb) containing SIAP3X was iso- 
lated from a male S. latifolia bacterial artificial chromosome (BAG) 
library, as described by Ishii et al. (2010), using the 4D-PGR method 
(Asakawa et al. 1997; Ishii et al. 2008). All BAG clones (13dllEa, 
13dllEb, and 7a8D) were purified using the PowerPrep HP Plasmid 
Purification Kit from OriGene Technologies (Rockville, MD) and 
were further purified by GsGl gradient ultracentrifugation. Emulsion 
PGR and 454 sequencing were performed on the Roche 454 FLX 
platform according to manufacturer's instructions. Sequence assembly 
was performed using the GS De novo Assembler. A summary of the 
sequencing data is shown in Supporting Information, Table SI. Gon- 
tigs consisting of more than 20 reads were first extracted and used to 
construct supercontigs in Sequencher (Gene Godes, Ann Arbor, MI). 
Gontigs of 13dllEa and 13dllEb were analyzed together. All gaps 
between the contigs were PGR-amplified and sequenced by the Sanger 
sequencing method. 

Sequence analysis 

Alignments of intron sequences of SIAP3X and SIAP3Y were con- 
ducted using the Needle program within the EMBOSS 6.5 package 
(Rice et al. 2000) with default settings, except that the gap penalty was 
set to 30.0 and the extend penalty was set to 0.0. Dot plot analyses 
were performed using HarrPlot 3.1.1 included in GENETYX version 
11 with the unit size to compare set to 10 and the dot plot matching 
number set to 8. Homology searches were performed with the nucle- 
otide sequence database of the National Genter for Biotechnology 
Information using the TBLASTX program and with the repetitive 
sequence database, Repbase (Jurka et al. 2005) of the Genetic 
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Information Research Institute (taxon: Viridiplantae), using the 
BLASTN program. These programs are contained in the BLAST+ 
package (Camacho et al 2009). Repetitive sequences, including LTRs, 
simple repeats, and other repeat sequences, were detected by homol- 
ogy searches of the 13dllE or 7a8D sequences with their own sequen- 
ces using the BLASTN program. The E-value threshold was set to 
IE— 10. Homologous regions were detected from homology searches 
of the 13dllE or 7a8D sequences using the BLASTN program. 
Thresholds for hit length, E-value, and identity were set to 500 bp, 
IE— 50, and 80%, respectively. The pair of coding regions of SIAP3X 
and SIAP3Y, the pairs of LTRs, and the pairs of homologous regions 
between the 13dllE and 7a8D sequences were aligned using the pro- 
gram MUSCLE (Edgar 2004). The synonymous nucleotide distance 
(k) between the coding regions of SIAP3X and S/AP3 7 was estimated 
using the Kumar method (Nei and Kumar 2000). The nucleotide 
distances (k) of the pairs of the LTRs and the pairs of the homologous 
regions were estimated using the Tamura-Nei model (Tamura and 
Nei 1993) within MEGA5 (Tamura et al 2011). An average substitu- 
tion rate (r) of 1.8 X 10~^ substitutions per synonymous site per year 
(Du et al 2010) was used for calculations. The time (T) since re- 
combination stopped between SIAP3X and SIAP3Y was estimated 
using the following formula: T = k/2r. 

Genomic Southern blot analysis 

Southern blot analysis methods were modified from those described 
previously (Kazama et al 2006). Genomic DNA was extracted from S. 
latifolia leaves using a Nucleon PhytoPure Genomic DNA Extraction 
Kit (GE Healthcare, Little Chalfont, England). Genomic DNA (15 |jLg) 
was digested with EcoRl for 12 hr. The concentrations of the digests 
were measured, and the digests were loaded onto a 1.0% (w/v) agarose 
gel and then transferred to an Immobilon-Ny-i- membrane (Merck, 
Darmstadt, Germany). Hybridization and signal detection were per- 
formed using the AlkPhos Direct Labeling and Detection System (GE 
Healthcare), and each probe was amplified individually by PGR with 
locus-specific primers (Table S2). The hybridized membranes were 
visualized using CDP-Stor (GE Healthcare) at room temperature with 
appropriate exposure times. 

FISH analysis 

Fluorescent in situ hybridization (FISH) was performed as previously 
described (Ishii et al 2010), with minor modifications. Two probes for 
the homologous (I) and nonhomologous (a) regions were prepared 
from the same PGR products produced in the genomic Southern blot 
analysis using DIG-Nick and the Biotin-Nick Translation Mix (Roche 
Diagnostics, Basel, Switzerland), respectively. Chromosomal DNA was 
denatured at 70° for 1 min in 2x SSC buffer containing 70% form- 
amide. The chromosome preparations were dehydrated immediately 
by 5-min treatments with 70% ethanol at —20° and 100% ethanol at 
room temperature. The preparations were then dried for 30 min at 
room temperature. Each slide was loaded with 10 |jlL of the hybrid- 
ization mixture containing 12.5 ng probe, 50% formamide, 10% dex- 
tran sulfate, and 2x SSC buffer. The slides were washed twice in 2x 
SSC buffer at 37°. The signals were fluoresceinated with antidigoxige- 
nin-rhodamine (Roche Diagnostics) and avidin-Alexa Fluor 488 (Mo- 
lecular Probes, Eugene, Oregon) for the digoxigenin-labeled probe and 
biotin-labeled probe, respectively. Preparations were counterstained 
with 4',6-diamidino-2-phenylindole (DAPI). The fluorescent Alexa 
Fluor 488, DAPI, and rhodamine images were detected consecutively 
under a Leica Q550 cytogenetic workstation (Leica Microsystems, 
Wetzlar, Germany) equipped with a black-and-white charged-coupled 
device camera (CoolSNAP HQ; Nippon Roper, Tokyo, Japan) using 



the SpectraVision filters SpectrumGreen, SpectrumDAPI, and Spec- 
trumOrange, respectively (Vysis, Chicago, IL). 

X chromosome linkage mapping 

To genotype the population, genomic PCRs were performed using K 
line-specific primers for SIAP3X, SlXl, and SIX4 (Table S3) on 96 F2 
progeny derived from a cross between a K line male and a B line 
female. Restriction fragment length polymorphism analysis was used 
to genotype the DD44X locus. Genomic DNA fragments were ampli- 
fied with DD44X-specific primers (Table S3) from the 96 F2 progeny, 
digested with Ncol, and analyzed by electrophoresis. The Kosambi 
mapping function was used to compute genetic distances in cM 
(Kosambi 1944). 

RESULTS 

Structure of the SIAP3X/Y transcriptional region 

The number of synonymous nucleotide divergences per synonymous 
sites between translated regions of SIAP3X and S/AP3 7 was 0.12 (SE, 
0.034). Nicolas et al (2005) estimated that in the oldest evolutionary 
strata, recombination stopped 5-10 MYA by using the average syn- 
onymous substitution rates of two genes (1.8 x 10 ~^ substitutions per 
synonymous site per year) in the Brassicaceae species (Koch et al 
2000). Similarly, we estimated that recombination between SIAP3X 
and SIAP3Y stopped 3.4 MYA (SE, 0.94). 

SIAP3X has five introns and SIAP3Y has six introns (Nishiyama 
et al 2010). Exon V and exon VI of SIAP3Y are homologous to the 
anterior half and the posterior half of exon V of SIAP3X (Nishiyama 
et al 2010), respectively, which indicated that intron 5 of S/AP37 was 
inserted after recombination stopped. Thus, the five pairs of introns 
shown in Table 1 were considered to be homologous because they 
underwent recombination before recombination suppression. Nish- 
iyama et al (2010) previously reported that intron 2 of SIAP3X was 
561 bp long, whereas intron 2 of SIAP3Y w^s 24.4 kb long, containing 
two retroelements and one telomere-like sequence. Therefore, we 
concluded that the majority of intron 2 of SIAP3Y was inserted after 
recombination stopped. We aligned each pair of introns and calcu- 
lated the similarities. Although lengths of the paired introns differed 
between SIAP3X and SlAP3Yy calculated sequence similarities (exclud- 
ing gaps) were between 76% and 82% (Table 1). 

We conducted dot plot analyses of each pair of introns (Figure 1). 
AU five pairs contained regions of sequence homology, which implied 
that the pairs were homologous before recombination suppression. 
However, these homologous regions did not span the full lengths of 
the introns but were interspersed like island chains. To investigate 
whether introns consisted of unique or repetitive sequences, we designed 
probes inside the introns and performed genomic Southern blot analyses 
(Figure 2). The probes designed for the introns showed three patterns: 
a strong smear pattern (introns 1 and 2 of SIAP3X and introns 1 and 4 of 
SIAP3Y); a weak smear pattern (intron 3, 4, and 5 of SIAP3X and introns 
3 of SIAP3Y); and a multiple-banded pattern (introns 5 and 6 of 
SIAP3Y). From this, it was revealed that all introns contained repetitive 
sequences; the distinctness of these patterns implied that each intron 
contained a different type of repeat sequence. 

Transposons in the peripheral regions of SIAP3X/Y 

To obtain sequences in the peripheral regions of SIAP3X and SIAP3Y, 
we sequenced the following three BAG clones using the 454-FLX 
sequencing platform: 13dllEa containing SIAP3X; 7a8D containing 
SIAP3Y (Nishiyama et al 2010); and 2b5E containing SIAP3X (iso- 
lated in this study). The sequences of 13dllEa and 2b5E spanned the 
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W Table 1 Homology between introns of SIAP3X and SIAP3Y 
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same chromosomal region, which indicated that they were derived 
from the same genome fragment; 2b5E was thus renamed 13dllEb. 
After the vector sequence was removed, we obtained a 11 5, 665 -bp 
sequence flanking SIAP3X (13dllE) and a 72,683-bp sequence flank- 
ing SIAP3Y (7a8D) (Figure 3). 

We conducted TBLASTX searches of the National Center for 
Biotechnology Information nucleotide database using the 13d HE and 
7a8D sequences as the query sequences (Table S4 and Table S5). In 
13d HE, 11 regions were homologous to the coding regions of LTR 
retrotransposons and one region was homologous to a coding region 
of a DNA transposon. In 7a8D, six regions were homologous to 
coding regions of LTR retrotransposons and two regions were ho- 
mologous to LINEs. BLASTN searches of the Genetic Information 
Research Institute Repbase database (taxon: Viridiplantae) using 
13dllE and 7a8D as query sequences revealed that 22 regions in 
13dllE were homologous to LTR retrotransposons, and seven 
regions in 7a8D were homologous to LTR retrotransposons (Table 
S6). The regions determined to be homologous to transposons by 
BLASTN were also identified using TBLASTX. Both 13dllE and 
7a8D sequences had regions homologous to gypsy retrotransposons 
and copia retrotransposons. Regions homologous to LINEs were found 
only in the 7a8D sequence. Genes other than transposons and SIAP3X/ 
y were not detected in either 13dllE or 7a8D. 

To identify LTR sequences sandwiching regions homologous to 
LTR retrotransposons, we performed pairwise BLASTN searches of 
13dllE and 7a8D within BACs (Figure 4). Six and four pairs of 



sequences, homologous in the forward direction, were found in 
13dllE and 7a8D, respectively, and were predicted to be LTR sequen- 
ces (Table 2 and Figure 4). All LTR retrotransposon-like elements did 
not contain any intact genes and were presumed to be nonautono- 
mous. At one and two loci in 13d HE and 7a8D, respectively, two LTR 
retrotransposon-like elements were identified with nested structures 
(Figure 4). The BLASTN searches also revealed the existence of two 
simple repeat sequences in 13dllE and seven simple repeat sequences 
in 7a8D. Six out of seven simple repeat sequences in 7a8D were found 
in one LTR retrotransposon-like element. The majority of sequences 
of the three LTR retrotransposon-like elements in 13dllE were ho- 
mologous to one another; similar results were observed for two DNA 
transposon-like elements in 13d HE, which indicated that they were 
genetically related to each other. 

Transposons inserted after recombination suppression 

Du et al. (2010) estimated the insertion dates of LTR retrotransposons 
using genetic distances calculated between left and right LTR sequen- 
ces. In the same manner, we calculated the genetic distances between 
the left and right LTR sequences of the LTR retrotransposon-like 
elements detected in 13dllE {SIAP3X) and 7a8D {SIAP3Y), from 
which we estimated the dates when they were inserted (Table 2 and 
Figure 5). AU the investigated LTR retrotransposon-like elements 
were estimated to be inserted after recombination between SIAP3X 
and SIAP3Y had stopped. Even the oldest one, which is the copia 
retrotransposon-like element, was presumed to be inserted 1.6 MYA 
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Figure 1 Dot plot analysis of 
introns of S\AP3X and S\AP3Y. 
Each intron sequence of S\AP3X 
(vertical) is compared to that of 
S/AP3y (horizontal). The first 700 
bp of S\AP3Y intron 2 is shown. 
Numbers above each axis indi- 
cate the ordinal number of each 
intron. Yellow squares indicate 
comparisons for which homolo- 
gous regions were detected. 
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Figure 2 Genomic Southern blot analysis of SIAP3X and SIAP3Y. (A) 
Schematics of genomic sequences of SIAP3X and SIAP3Y. White and 
black boxes indicate untranslated regions and exons, respectively. 
Numbers indicate intron numbers of the genes and bold lines indicate 
probe sequence coverage of each intron for Southern hybridization. 
(B) Genomic Southern blot analysis of introns of SIAP3X and SIAP3Y. 
Numbers above lanes indicate intron numbers of the genes corre- 
sponding to those in (A). 

(SE, 0.50). The dates when the LTR retrotransposon-like elements 
were inserted differed depending on the kind of element. The two 
copia retrotransposon-like elements in the 13dllE sequence were 
assumed to be inserted 1.2 MYA (SE, 0.48) and 1.6 MYA (SE, 
0.50). The one and two gypsy retrotransposon-like elements in the 
13dllE and 7a8D sequences were thought to be inserted 0.22 MYA 
(SE, 0.16), 0.62 MYA (SE, 0.11), and 0.65 MYA (SE 0.12), respectively. 
The three unknown retrotransposon-like elements in 13dllE were 
considered to be inserted 0.18 MYA (SE, 0.057), 0.30 MYA (SE, 
0.079), and 0.74 MYA (SE 0.12), respectively. One and two pairs of 
nested LTR retrotransposon-like elements were detected in 13dllE 
and 7a8D, respectively (Figure 5). It was natural to assume that 
the inner element was inserted after the outer element. With regard 
to the two pairs in the 7a8D sequence, the estimated insertion dates of 
the inner elements were more recent than those of the outer elements. 



The estimated insertion time of the inner element of the pair identified 
in 13dllE was older than that of the outer element, but the two in- 
sertion dates were determined to be close and not significantly different. 

Homology between noncoding DNA regions of the X 
and Y chromosomes 

The noncoding DNA regions around SIAP3X and S/AP3 l^except for 
the transposons that were inserted after recombination stopped — 
should have been homologous. We attempted to reveal structural 
changes of the sex chromosomes after recombination suppression 
by investigating how much homologous sequence remained. To detect 
homologous regions between the X and Y chromosomes, we con- 
ducted BLASTN searches of 13dllE (SIAP3X) using 7a8D (SIAP3Y) 
as the query sequence. Homologous regions were scattered around 
SIAP3X/Y (Table 3 and Figure 6A). The total length of the homolo- 
gous regions in the 13dllE and 7a8D sequences were 7.7 kb (6.7%) 
and 7.7 kb (10.6%), respectively. We calculated the nucleotide diver- 
gence of each homologous sequence pair (Table 3). The divergences 
showed some variation, ranging between 0.085 (SE, 0.0085) (Table 3, 
ID 4) and 0.25 (SE, 0.029) (Table 3, ID 5), whereas that observed for 
the coding regions of SIAP3X/Y w^s 0.12 (SE, 0.034). 

To investigate whether the homologous regions between 13dllE 
and 7a8D in the noncoding DNA regions were unique in the genome, 
we performed genomic Southern hybridizations with probes made 
from each PGR product of the homologous region in the 13dllE 
sequence (Table 3 and Figure 6B). Although probe sequences do not 
include transposons, all the hybridizations produced smeary patterns, 
which implied that these sequences were abundant in the genome. 
FISH performed using one of those sequences (Figure 6B, I) as a probe 
also yielded many signals on multiple chromosomes (Figure 6C). To 
investigate the uniqueness of noncoding DNA regions in the genome, 
which did not show homology to their counterparts on the opposite 
chromosome (Figure 6 A, horizontal lines), we also analyzed genomic 
Southern hybridization data using PGR products generated from these 
DNA regions as probes (Figure 6B, a-j). All hybridizations yielded 
smeared patterns, which indicated that the genome included a large 
number of these sequences. FISH performed using one of these 
sequences (Figure 6B, d) as a probe also showed signals on multiple 
chromosomes, which supported the repetitive nature of the sequence 
(Figure 6D). 

DISCUSSION 

Beginning of divergence between sequences 
surrounding SIAP3X and SIAP3Y 

We estimated that recombination between SIAP3X and SIAP3Y stop- 
ped 3.4 MYA. In addition, we determined the dates of the subsequent 
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Figure 3 Schematics of BAC 
sequences. Horizontal lines in- 
dicate entire sequences of each 
BAC clone drawn to scale as 
indicated. Black bars indicate 
exons, which are labeled with 
Roman numerals. 
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Figure 4 Transposable elements in the BAC sequences. Entire sequences of the BAC clones and the locations of exons are displayed in the same 
manner as in Figure 3. Blue, green, gray, and orange boxes on the horizontal line indicate regions homologous to copia-like retrotransposons, 
gypsy-like retrotransposons, DNA transposons, and LINEs, respectively. Bold and thin boxes are detected by the TBLASTX and BLASTN searches, 
respectively. White boxes sandwiched by pairs of red boxes connected by arcs on the horizontal line indicate the locations of LTR retrotransposon-like 
elements. The numbers below LTRs correspond to those listed in Table 2. Brown boxes indicate simple repeat sequences. Pairs of boxes outlined in 
solid red lines and connected by an arc above the horizontal line indicate other repetitive sequences. Pairs of boxes outlined with dashed red lines 
indicate homologous regions in which evolutionarily related elements have been inserted. 



insertions of transposons in this region. The molecular clock used for 
age estimations, which was the same as the one Nicolas et al. (2005) 
used, was based on the Brassicaceae species, which is currently the 
most reliable clock for S. latifolia. More accurate estimation will be 
possible when a molecular clock of more closely related species 
becomes available. 

To confirm the location of SIAP3X within the nonrecombining 
region, we constructed a linkage map of the S. latifolia X chro- 
mosome based on the estimation of genetic distances (Table S7). 
The most probable order of the four genes from the PAR on the X 
chromosome was S/Xi, SIAP3X, DD44X, and SIX4. The relative 
locations of genes on the X chromosome, integrated with the data 
presented by Bergero et al. (2013), are shown in Figure SI. Twenty 
percent divergence is typical for stratum 1 genes and 10% diver- 
gence is typical for stratum 2 genes (Cegan et al 2010). SIAP3X 
was mapped between DD44X and SlXly and the divergence be- 
tween SIAP3X and Y was 12%, which implies that SIAP3X is lo- 
cated in stratum 2. 



Differentiation of sequences in sex chromosomes 

Many diverged transposons were found in 13dllE {SIAP3X) and in 
7a8D (SIAP3Y) (Figure 4); however, the sequence contents of these 
elements were determined to be different. For example, no significant 
similarity was found between retroelements identified in 13dllE 
iSlAP3X) and 7a8D {SIAP3Y), indicating either that these retroele- 
ments were independently inserted into 13dllE {SIAP3X) or 7a8D 
{SIAP3Y) after recombination stopped or that they had existed on the 
proto-sex chromosomes and were translocated to a distant region on 
only one of the chromosomes. Because the estimated insertion time of 
the oldest LTR retroelement (1.6 MYA) is more recent than the pre- 
dicted time when recombination between SIAP3X and SIAP3Y (3.4 
MYA) stopped, the former scenario is more likely to have occurred. 
Gschwend et al (2012) compared the sequence of the C. papaya X 
chromosome with the autosomal counterpart of the related species, 
Vasconcellea monoica, and reported that sequences not only on the Y 
chromosome but also on the X chromosome expanded after recom- 
bination suppression. Obtaining the sequence of the SIAP3 ortholog 



W Table 2 Predicted LTR retroelements and estimated insertion times 



Clone 


ID 


Type 


Start, bp 


End, bp 


Length, bp 


Left LTR 
Length, bp 


BLASTN 
E-value Identity, % 


K(SE) 


Estimated Insertion 
Time, MYA (SE) 


13d11 E (S/AP3X) 


1 


Unknown 


8721 


14,676 


5956 


1548 


0 


99 


0.0065 (0.0021) 


0.18 (0.057) 




2 


Unknown 


23,393 


29,638 


6246 


1387 


0 


99 


0.011 (0.0028) 


0.30 (0.079) 




3 


Gypsy 


54,004 


70,569 


10,886 


715 


0 


99 


0.0080 (0.0057) 


0.22 (0.16) 




4 


Unknown 


58,455 


64,134 


5680 


1450 


0 


96 


0.027 (0.0044) 


0.74 (0.12) 




5 


Copia 


74,749 


79,514 


4766 


169 


7E-65 


93 


0.045 (0.017) 


1 .2 (0.48) 




6 


Copia 


88,141 


92,782 


4642 


208 


5E-86 


95 


0.057 (0.018) 


1.6 (0.50) 


7a8D {SIAP3y) 


7 


Unknown 


12,651 


32,108 


14,093 


366 


4E-153 


94 


0.044 (0.012) 


1.2 (0.32) 




8 


Unknown 


21,888 


27,252 


5365 


1670 


0 


99 


0.0030 (0.0013) 


0.083 (0.037) 




9 


Gypsy 


36,730 


58,367 


12,663 


1538 


0 


97 


0.022 (0.0039) 


0.62 (0.11) 




10 


Gypsy 


45,305 


54,279 


7770 


1345 


0 


97 


0.024 (0.0043) 


0.65 (0.12) 



K, nucleotide divergence between left and right LTRs; BLASTN, homology between the left and right LTR sequences calculated by the BLASTN program. 



2126 I K.lshii eta/. jli: G3- Genes | Genomes | Genetics 



SIAP3X {^36^^E) 
0 r X ■ X3I ■ El ^ 




SIAP3Y {7aQD) 



Figure 5 Estimated insertion 
dates of LTR retroelements. 
The horizontal line in the middle 
of the vertical axis indicates the 
constructed genomic sequence 
peripheral to the ancient SIAP3X/ 
Y before recombination suppres- 
sion. The upper half and lower 
half of the figure show temporal 
changes of genomic sequence 
peripheral to S/AP3X and SIAP3Y, 
respectively, by insertion of LTR 
retroelements in single logarith- 
mic plots. LTR retroelements are 
displayed in the same manner as 
in Figure 4. Horizontal lines drawn 
through the center of boxes indi- 
cate the estimated insertion date 
of LTR retroelements. Dashed 
lines indicate insertion sites. 



and its flanking region in Silene vulgaris, which is the closest related 
hermaphroditic species to S. latifolia, would enable us to better discern 
changes in the DNA sequence that occurred independently in the X 
and Y chromosomes. 

AU three unknown retrotransposon-like elements (ID 1, ID 2, and 
ID 4 in Table 2 and Figure 4) in 13dllE {SIAP3X) showed evidence of 
homology. In addition, the sequences of the two DNA transposon-like 
elements, which were nested in 13dllE (S/AP3X), were also homolo- 
gous to one another (Figure 4). Thus, each group of transposons was 
concluded to have the same origin. The estimated insertion times of 
the three unknown retrotransposon-like elements in 13dllE {SIAP3X) 
were similar (0.18-0.74 MYA) (Table 2), suggesting that transposons 
of this type were translocated wdthin a limited period. The estimated 
insertion times of the two copia retrotransposon-like elements in 
13dllE {SIAP3X) were also similar (1.2 and 1.6 MYA) (Table 2). It 
is possible that each group of transposons has a shared origin but 
has lost homology over time. LINEs were detected only in the 7a8D 
sequence (Figure 4); however, according to the mapping of trans- 
posons by FISH analyses (Kejnovsky et al. 2009), they were distrib- 
uted across the full length of both the X and Y chromosomes. It is 
possible that LINEs could be identified in the adjacent sequence of 
13dllE {SIAP3X). 

Yu et al (2008) compared two pairs of BAG sequences derived 
from sex chromosomes of C. papaya. One pair of sequences (X BAG 
61H02 and Y^ BAG 95B12) was thought to have stopped recombining 
7.0 MYA, in that 81.6% and 74.5% of sequences of the X and Y^ 



chromosome had homology, respectively, and the identity between 
homologous sequences was 87.4%. The other pair of sequences (X 
BAG 53E18 and Y^ BAG 85B24) was estimated to have ceased re- 
combination 1.9 MYA. In this example, 41.5% and 30.7% of sequences 
of the X and Y^ chromosome had homology, respectively, and the 
identity between homologous sequences was 83.6%. In the pair of 
sequences (13dllE and 7a8D) analyzed here for S. latifolia, which 
were assumed to have stopped recombining 3.4 MYA, 6.7% and 
10.6% of sequences of the X and Y chromosomes had homology, 
respectively, and the average identity between homologous sequen- 
ces was 76.9%. Differences in the percentages of the noncoding 
homologous region are thought to be, in part, attributable to differ- 
ences in the lengths of the analyzed sequences. The lower percen- 
tages despite the identity may reflect a high rate of accumulation of 
transposable elements, which may have partially caused the hetero- 
morphism of sex chromosomes in this species. 

Regarding indels, Sundstrom et al (2003) conducted comparative 
analysis on intron sequences of paralogs in primates and birds and 
revealed that in birds with female heterogamety, indels occurred in the 
Z chromosome approximately twice as often as in the W chromo- 
some, and that in primates with male heterogamety, indels occurred at 
similar rates in both the X and Y chromosome. This indicated that 
both meiotic recombination and the number of cell divisions in re- 
productive cells play an important role in the occurrence of indels. It 
raises an interesting question regarding whether this theory is also 
applicable to plants in which reproductive cells are generated from 



■ Table 3 Homologous regions between peripheral regions of SIAP3X and SIAP3Y 

1 3d11 E {SIAP3X) 7a8D {SIAP3Y) 



ID 


Start, bp 


End, bp 


Length, bp 


Start, bp 


End, bp 


Length, bp 


Direction 


E-value 


Identity, % 


K(SE) 


1 


3586 


5587 


2002 


71,985 


70,005 


1981 


Inverted 


0 


73 


0.14 (0.012) 


2 


16,253 


17,744 


1492 


8055 


9530 


1474 


Direct 


0 


76 


0.11 (0.010) 


3 


17,735 


18,992 


1258 


3248 


1996 


1253 


Inverted 


0 


81 


0.14 (0.014) 


4 


31,730 


33,237 


1508 


1996 


3467 


1470 


Direct 


0 


78 


0.085 (0.0085) 


5 


96,491 


97,215 


725 


67,473 


66,715 


759 


Inverted 


5.00E-166 


78 


0.25 (0.029) 


6 


96,492 


97,219 


728 


13,016 


13,789 


772 


Direct 


4.00E-161 


77 


0.092 (0.013) 



K, nucleotide divergence between the regions in 13d11E and 7a8D. 
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Figure 6 Noncoding DNA regions in peripheral regions of SIAP3X and SIAP3Y. (A) Maps of noncoding DNA regions. Entire sequences of the 
BAC clones and the locations of exons are displayed in the same manner as in Figure 3. White boxes indicate coding regions of transposable 
elements or regions sandwiched by a pair of LTRs. Pink and blue boxes indicate noncoding and homologous regions between 1 3d1 1 E and 7a8D 
in the forward direction and reverse direction when S/AP3X and S/AP3y are aligned in the same direction, respectively. Numbers on the boxes 
indicate correspondence between homologous regions. Horizontal lines indicate regions that are noncoding but not homologous between 13d11E 
and 7a8D. Black bold lines indicate probe regions for genomic Southern hybridization and FISH. Roman numbers and lower-case characters below 
bold lines correspond to those in (B), (C), and (D). (B) Genomic Southern blot analysis of noncoding DNA regions. The membranes used in Figure 2B 
were hybridized with each PGR product of noncoding DNA regions. Roman and Arabic numbers above lanes indicate probe regions corresponding 
to those in (A). (G) FISH result for the probe of the region (I) that covered the noncoding and X-Y homologous sequence. (D) FISH result for the probe 
of the region (a) spanning the noncoding and nonhomologous sequence. 



vegetative cells, thus weakening the influence of the difference in the 
number of cell divisions between male and female reproductive cells. 
However, regarding transposable elements, a comparison between the 
40-kb neo-X and 45-kb neo-Y sequences of Drosophila miranda, 
which is assumed to have emerged approximately 1 MYA, revealed 
that insertions of seven transposable elements occurred only in the 
neo-Y chromosome (Bachtrog 2005). It is assumed that transposable 
elements tend to accumulate more in the recombination suppression 
region of the Y chromosome. In the regions in S. latifolia analyzed in 
this study, transposable elements have highly accumulated not only in 
the Y chromosome but also in the X chromosome for 3.4 million 
years, preventing us from revealing whether the insertion of the trans- 
posable element was Y chromosome-biased. Obtaining continuous 
sequences between SIAP3 and its neighboring genes in both X and 
Y chromosomes will enable us to find the answers to our questions. 

Origin of the noncoding DNA region 

AU sequences of noncoding DNA regions investigated in this study 
contained repetitive elements, irrespective of whether they were homol- 
ogous between sex chromosomes (Figure 2B and Figure 6, B-D). Five 
noncoding DNA regions exhibited homology between the X and Y 
chromosomes (Table 3, IDs 1-4 and ID 6) and showed no significant 



differences in divergence from the synonymous nucleotide divergences 
observed for the translated regions of SIAP3X and SIAP3Y, as determined 
by the chi-square test. Because the divergence estimates for these regions 
along with SIAP3X/Y are consistent with their chromosomal location in 
stratum 2, it is considered that these regions are derived from ancestral 
chromosomes. However, one noncoding DNA region (Table 3, ID 5) had 
higher divergences (P = 0.0001 by chi-square test). Although the cause of 
the high divergence is uncertain, it is presumed that the sequences of the 
region have paralogous origins that had been duplicated in the past. 
Alternatively, as suggested in primates (Iwase et al 2010), it could be 
possible that gene conversion occurred locally between the SIAP3X and 
SIAP3Y coding regions after recombination suppression, which would 
make divergence lower. With regard to noncodiag DNA regions that 
were not homologous between the sex chromosomes, it was assumed 
that these sequences were inherited from ancestral chromosomes and 
that their counterparts in either the X chromosome or the Y chromo- 
some were lost during deletion events. 

Evolutionary strata on the X chromosome in 5. latifolia 

Bergero et al (2013) revealed that the X chromosome of S. latifolia 
have three evolutionary strata. However, because whole sequences of 
sex chromosomes have not been generated, it is uncertain whether the 
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evolutionary strata on the X chromosome of S. latifolia have clear 
boundaries. Lemaitre et al (2009) proposed that evolutionary strata 
on the human X chromosome were formed by inversions on the Y 
chromosome, which has been explained by the isochromatid model 
(Ranz et al 2007), because boundary sequences between strata 4 and 
strata 5 and between PAR and strata 5 on the X chromosome are 
duplicated on the Y chromosome. According to the isochromatid 
model, two pairs of staggered single- strand breaks yield long 5' over- 
hangs, which are filled-in by DNA synthesis. When such breaks are 
repaired by nonhomologous end-joining, they can result in inversions 
flanked by inverted duplications of the sequences between the paired 
single-strand breaks (Ranz et al 2007). 

In two instances, we identified noncoding DNA regions that were 
present at two locations on one sex chromosome and only one 
homologous location on the second sex chromosome (Figure 6A). For 
example, region 3/4 in 7a8D {SIAP3Y) was homologous to region 3 
and region 4 in 13dllE {SIAP3X). Region 2 in 7a8D, which was 
adjacent to region 3/4, was also homologous to region 2 in 13dllE, 
which was adjacent to region 3 (Figure 6A). These homologous non- 
coding DNA regions could be vestiges of inversions; however, defin- 
itive boundaries of sequence similarity that could be used to indicate 
evolutionary strata were not found. 

In this study, we analyzed 116 kb and 73 kb of sequence from the X 
and Y chromosomes of S. latifolia and found that only 6.7% of the X 
chromosome-derived sequence and 10.6% of the Y chromosome-de- 
rived sequence were homologous between the two chromosomes. It was 
shown that sex chromosomes of S. latifolia were heteromorphic not only 
at the chromosomal level but also at the sequence level. We detected two 
locations in which homologous sequences between the X and Y chro- 
mosomes were duplicated, implying the occurrence of inversions. A 
wider comparison of noncoding DNA regions between the X and Y 
chromosomes, especially those near the PAR, will help delineate the 
boundaries of evolutionary strata and define differences in the degree 
of heteromorphy between these strata, which will be essential for gaining 
a broader understanding of the evolution of plant sex chromosomes. 

ACKNOWLEDGMENTS 

We thank Yoko Kuroki and Asao Fujiyama for technical support as 
a part of the support research project conducted by MEXT KAKENHI 
"Genome Science" (grant 221S0002). This work was supported by 
JSPS KAKENHI (grant 09J07377 to K.I., grants 20780009, 23770070, 
and 25292009 to Y.K., and grants 18657001 and 24657046 to S.K.), by 
incentive research grants from RIKEN (to K.I.), by the Fund for Seeds 
of Collaborative Research from RIKEN (to Y.K.). 

LITERATURE CITED 

Asakawa, S., I. Abe, Y. Kudoh, N. Kishi, Y. Wang et al, 1997 Human BAG 
library: construction and rapid screening. Gene 191: 69-79. 

Bachtrog, D., 2005 Sex chromosome evolution: Molecular aspects of Y- 
chromosome degeneration in Drosophila. Genome Res. 15: 1393-1401. 

Bergero, R., A. Forrest, E. Kamau, and D. Gharlesworth, 2007 Evolutionary 
strata on the X chromosomes of the dioecious plant Silene latifolia: ev- 
idence from new sex-linked genes. Genetics 175: 1945-1954. 

Bergero, R., S. Qiu, A. Forrest, H. Borthwick, and D. Gharlesworth, 

2013 Expansion of the pseudo-autosomal region and ongoing recom- 
bination suppression in the Silene latifolia sex chromosomes. Genetics 
194: 673-686. 

Gamacho, G., G. Goulouris, V. Avagyan, N. Ma, J. Papadopoulos et al, 

2009 BLAST-h: architecture and applications. BMG Bioinformatics 10: 421. 
Gegan, R., G. A. B. Marais, H. Kubekova, N. Blavet, A. Widmer et al, 

2010 Structure and evolution of ApetalaS, a sex-linked gene in Silene 
latifolia. BMG Plant Biol. 10: 180. 



Gharlesworth, D., 2002 Plant sex determination and sex chromosomes. 

Heredity 88: 94-101. 
Gharlesworth, D., 2013 Plant sex chromosome evolution. J. Exp. Bot. 64: 

405-420. 

Gharlesworth, D., and B. Gharlesworth, 1980 Sex differences in fitness and 
selection for centric fusions between sex-chromosomes and autosomes. 
Genet. Res. 35: 205-214. 

Du, J., Z. Tian, G. S. Hans, H. M. Laten, S. B. Gannon et al, 

2010 Evolutionary conservation, diversity and specificity of LTR-ret- 
rotransposons in flowering plants: insights from genome-wide analysis 
and multi-specific comparison. Plant J. 63: 584-598. 

Edgar, R. G., 2004 MUSGLE: multiple sequence alignment with high ac- 
curacy and high throughput. Nucleic Acids Res. 32: 1792-1797. 

Gschwend, A. R, Q. Yu, E. J. Tong, F. Zeng, J. Han et al, 2012 Rapid 
divergence and expansion of the X chromosome in papaya. Proc. Natl. 
Acad. Sci. USA 109: 13716-13721. 

Ishii, K., R. Sugiyama, M. Onuki, Y. Kazama, S. Matsunaga et al, 2008 The 
Y chromosome-specific STS marker MS2 and its peripheral regions on 
the Y chromosome of the dioecious plant Silene latifolia. Genome 51: 
251-260. 

Ishii, K., Y. Amanai, Y. Kazama, M. Ikeda, H. Kamada et al, 2010 Analysis 
of BAG clones containing homologous sequences on the end of the Xq 
arm and on chromosome 7 in the dioecious plant Silene latifolia. Genome 
53: 311-320. 

Iwase, M., Y. Satta, Y. Hirai, H. Hirai, H. Imai et al, 2003 The amelogenin 
loci span an ancient pseudoautosomal boundary in diverse mammalian 
species. Proc. Natl. Acad. Sci. USA 100: 5258-5263. 

Iwase, M., Y. Satta, H. Hirai, Y. Hirai, and N. Takahata, 2010 Frequent gene 
conversion events between the X and Y homologous chromosomal re- 
gions in primates. BMG Evol. Biol. 10: 225. 

Jurka, }., V. V. Kapitonov, A. Pavlicek, P. Klonowski, O. Kohany et al, 
2005 Repbase Update, a database of eukaryotic repetitive elements. 
Gytogenet. Genome Res. 110: 462-467. 

Kazama, Y., R. Sugiyama, S. Matsunaga, F. Shibata, W. Uchida et al, 

2003 Organization of the Kpnl family of chromosomal distal-end sat- 
ellite DNAs in Silene latifolia. J. Plant Res. 116: 317-326. 

Kazama, Y., R. Sugiyama, Y. Suto, W. Uchida, and S. Kawano, 2006 The 
clustering of four subfamilies of satellite DNA at individual chromosome 
ends in Silene latifolia. Genome 49: 520-530. 

Kejnovsky, E., R. Hobza, T. Germak, Z. Kubat, and B. Vyskot, 2009 The role 
of repetitive DNA in structure and evolution of sex chromosomes in 
plants. Heredity 102: 533-541. 

Koch, M., B. Haubold, and T. Mitchell-Olds, 2000 Gomparative evolu- 
tionary analysis of chalcone synthase and alcohol dehydrogenase loci in 
Arahidopsis, Arabis and related genera (Brassicaceae). Mol. Biol. Evol. 17: 
1483-1498. 

Kosambi, D. D., 1944 The estimation of map distance from recombination 
values. Ann. Eugen. 12: 172-175. 

Lahn, B. T., and D. G. Page, 1999 Four evolutionary strata on the human X 
chromosome. Science 286: 964-967. 

Lawson-Handley, L. }., H. Geplitis, and H. Ellegren, 2004 Evolutionary 
strata on the chicken Z chromosome: implications for sex chromosome 
evolution. Genetics 167: 367-376. 

Lemaitre, C, M. D. Braga, G. Gautier, M. F. Sagot, E. Tannier et al, 
2009 Footprints of inversions at present and past pseudoautoso- 
mal boundaries in human sex chromosomes. Genome Biol. Evol. 1: 
56-66. 

Liu, Z., P. H. Moore, H. Ma, G. M. Ackerman, M. Ragiba et al, 2004 A 
primitive Y chromosome in papaya marks incipient sex chromosome 
evolution. Nature 427: 348-352. 

Matsunaga, S., M. Hizume, S. Kawano, and T. Kuroiwa, 1994 Gytological 
analyses in Melandrium album: genome size, chromosome size and 
fluorescence in situ hybridization. Gytologia (Tokyo) 59: 135-141. 

Matsunaga, S., S. Kawano, H. Takano, H. Uchida, A. Sakai et al, 

1996 Isolation and developmental expression of male reproductive 
organ-specific genes in a dioecious campion, Melandrium album {Silene 
latifolia). Plant J. 10: 679-689. 



r-^: GS'Genes | Genomes | Genetics 



Volumes December 2013 I Noncoding DNA in S. latifolia Sex Chromosomes I 2129 



Matsunaga, S., S. Kawano, and T. Kuroiwa, 1997 MROSl, a male stamen- 
specific gene in the dioecious campion Silene latifolia is expressed in 
mature pollen. Plant Cell Physiol. 38: 499-502. 

Matsunaga, S., E. Isono, E. Kejnovsky, B. Vyskot, J. Dolezel et al, 

2003 Duplicative transfer of a MADS box gene to a plant Y chromo- 
some. Mol. Biol. Evol. 20: 1062-1069. 

Ming, R., S. Hou, Y. Feng, Q. Yu, A. Dionne-Laporte et al, 2008 The draft 
genome of the transgenic tropical fruit tree papaya {Carica papaya Lin- 
naeus). Nature 452: 991-996. 

Mosiolek, M., P. Pasierbek, J. Malarz, M. Mos, and A. }. Joachimiak, 
2005 Rumex acetosa Y chromosomes: constitutive or facultative het- 
erochromatin? Folia Histochem. Cytobiol. 43: 161-167. 

Nam, K., and H. EUegren, 2008 The chicken {Gallus gallus) Z chromosome 
contains at least three non-linear evolutionary strata. Genetics 180: 1131-1136. 

Nei, M., and S. Kumar, 2000 Molecular Evolution and Phylogenetics, Oxford 
University Press, Oxford, England. 

Nicolas, M., G. Marais, V. Hykelova, B. Janousek, V. Laporte et al, 2005 A 
gradual process of recombination restriction in the evolutionary history 
of the sex chromosomes in dioecious plants. PLoS Biol. 3: e4. 

Nishiyama, R., K. Ishii, E. Kifune, Y. Kazama, K. Nishihara et al, 2010 Sex 
chromosome evolution revealed by physical mapping of SIAP3X/Y in the 
dioecious plant Silene latifolia. Cytologia (Tokyo) 75: 319-325. 

Ranz, J. M., D. Maurin, Y. S. Chan, M. von Grotthuss, L. W. Hillier et al, 
2007 Principles of genome evolution in the Drosophila melanogaster 
species group. PLoS Biol. 5: el 52. 

Renner, S. S., and R. E. Ricklefs, 1995 Dioecy and its correlates in the 
flowering plants. Am. J. Bot. 82: 596-606. 

Rice, P., I. Longden, and A. Bleasby, 2000 EMBOSS: The European Mo- 
lecular Biology Open Software Suite. Trends Genet. 16: 276-277. 



Siroky, J., M. A. Lysak, J. Dolezel, E. Kejnovsky, B. Vyskot et al, 

2001 Heterogeneity of rDNA distribution and genome size in Silene 
spp. Chromosome Res. 9: 387-393. 

Sundstrom, H., M. T. Matthew, and H. EUegren, 2003 Insertion and de- 
letion mutation male biased?: Molecular evolutionary analysis of avian 
and primate sex chromosome sequences. Genetics 164: 259-268. 

Tamura, K., and M. Nei, 1993 Estimation of the number of nucleotide 
substitutions in the control region of mitochondrial DNA in humans and 
chimpanzees. Mol. Biol. Evol. 10: 512-526. 

Tamura, K., D. Peterson, N. Peterson, G. Stecher, M. Nei et al, 

2011 MEGA5: Molecular evolutionary genetics analysis using maxi- 
mum likelihood, evolutionary distance, and maximum parsimony 
methods. Mol. Biol. Evol. 28: 2731-2739. 

Telgmann-Rauber, A., A. Jamsari, M. S. Kinney, }. C. Pries, and C. Jung, 
2007 Genetic and physical maps around the sex-determining M- 
locus of the dioecious plant asparagus. Mol. Genet. Genomics 278: 
221-234. 

Wang, J., J. K. Na, Q. Yu, A. R. Gschwend, J. Han et al, 2012 Sequencing 
papaya X and Y^ chromosomes reveals molecular basis of incipient 
sex chromosome evolution. Proc. Natl. Acad. Sci. USA 109: 13710- 
13715. 

Yamato, K. T., K. Ishizaki, M. Fujisawa, S. Okada, S. Nakayama et al, 

2007 Gene organization of the liverwort Y chromosome reveals distinct 
sex chromosome evolution in a haploid system. Proc. Natl. Acad. Sci. 
USA 104: 6472-6477. 

Yu, Q., S. Hou, F. A. Feltus, M. R. Jones, J. E. Murray et al, 2008 Low X/Y 
divergence in four pairs of papaya sex-linked genes. Plant J. 53: 124-132. 

Communicating editor: N. D. Young 



2130 I K. Ishii eta/. 



r-^.G3' Genes \ Genomes | Genetics 



